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Abstract 



A study of rotational properties of the ground super deformed bands in 190 Hg, 
192 Hg, 194 Hg, and 194 Pb is presented. We use the cranked Hartree-Fock-Bogoliubov 
method with the SkM* parametrization of the Skyrme force in the particle-hole 
channel and a seniority interaction in the pairing channel. An approximate parti- 
cle number projection is performed by means of the Lipkin-Nogami prescription. 
We analyze the proton and neutron quasiparticle routhians in connection with the 
present information on about thirty presently observed superdeformed bands in 
nuclei close neighbours of 192 Hg. 



1 Introduction 

With the recent evolution in 7-ray multidetectors, the experimental knowledge of super- 
deformed (SD) bands in nuclei of the mercury region is rapidly growing. Since the initial 
discovery of a SD rotational band in 191 Hg |lj, more than forty bands have been observed 
for isotopes of Au, Hg, Tl and Pb (see refs. §, g 1, g |, 0, |, |, 0, HUl and the earlier 



papers cited in the compilation of experimental data The anchor point of the region 

appears to be the nucleus 192 Hg which corresponds to the SD magic numbers: Z=80 and 
A=112. 

Although superdeformation has been predicted by many models (for a review see 
ref. [0), the spectrum of approaches used up to now to study high spin properties of 
the SD states is still rather limited. The theoretical analyses have been mostly done in 
terms of the cranked Nilsson-Strutinsky Jl5[ method, either with the modified harmonic 



oscillator (MHONS), e.g. ref. pj, or with the Wood-Saxon (WSNS), e.g. ref. |]| |l| 
potentials. The methods relying on a mean field derived from a microscopic nucleon- 
nucleon force have been applied only in two cases. For 192 Hg, 194 Hg and 194 Pb the cranked 
Hartree-Fock equations have been solved [2(J and in 192 Hg the Bohr hamiltonian approach 



with the collective potential and mass parameters calculated at zero spin has been used 
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In view of their general good quality at describing zero-spin properties, the latter 
methods should be more extensively tested at high spins. Such methods are free from 
geometric assumptions concerning a detailed parametrization of the nuclear surface, and 
therefore they seem to be more apt at describing polarization phenomena that may appear 
in the dynamical behaviour of a rotating nucleus. 

As compared to the other major region of superdeformation (around v4=150) the 
high spin properties of nuclei in the 192 Hg region seem to be different because of the 
presence of large pairing correlations. In particular, the generally observed increase of 
the moment of inertia is interpreted |22|, |23| as resulting both from the alignment of 
intruder orbitals and from a gradual disappearance of pairing correlations. The relative 
importance of these two effects requires further microscopic investigations. There exist 
also uncertainties concerning the predictions of crossing frequencies, the determination 
of pairing strengths and the evaluation of the relative importance of dynamical pairing 
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correlations compared to those described by the static pairing formalism which is usually 
sufficient to account for zero spin properties. 

The present work is devoted to a presentation of the cranked Hartree-Fock-Bogoliubov 
(HFB) method and to a discussion of its ability to address such issues. The next section 
contains a short description of our theoretical approach. In particular, we present a 
method of solution of the cranked HFB equations with the Skyrme interaction in the 
coordinate representation. Pairing correlations are described by a seniority interaction 
and we briefly describe how we implement an approximate variation after projection by 
means of the Lipkin-Nogami approach |24|, |25], ^6f (HFBLN). 

As a first application we use the cranked HFB method to study the SD bands in 
nuclei around 192 Hg. Section 3 contains a presentation of our results. First we investigate 
properties of the ground SD bands in 190 Hg, 192 Hg, 194 Hg and 194 Pb. Then, we show 
that the properties of SD bands observed in other odd and even nuclei in this region 
are compatible with our spectra of quasiparticle routhians. The last section presents our 
conclusions. 



2 The HFBLN method 



2.1 THE TWO-BASIS METHOD 



First, we review the cranked HFB method [27] as much as required to introduce our 
notation. Given the hamiltonian H which in a single particle basis {af, a,i,(i = l,..., N)} 
has the form 

H = f+V , 

T = 'Yliij T^a i cij , (1) 

V = 4 J2ijkl VijklO>t a t a l a k j 

this method determines the quasiparticle vacuum state which minimizes the energy 

£ = (V\H\9) (2) 

subject to the constraints (\l/|iV|\I/) = N and J x \^) = Jo- The average value of 
the particle-number operator N is thus constrained to the desired number of nucleons 
iVo (separately for neutrons and protons) and that of the cranking operator J x to the 
desired value of angular momentum^ Jo. The collective-cranking operator J x is equal 
to the projection of the angular-momentum vector J on the axis perpendicular to the 
largest-elongation axis. 

The determination of amounts to finding the unitary matrix W which relates 
the particle creation and annihilation operators {af, <2j} to the quasiparticle operators 

1 In general, Jo is related to the angular momentum J and its projection K on the symmetry axis by 
Jo = \J J( J + 1) — K 2 . Since we are mostly interested in high spin properties and K—0 bands we take 
J «J. 
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{bf, bi} such that the fo's are destruction operators associated with the quasiparticle 
vacuum |*) (foil*) = 0, % = 1, . . . , N): 

The minimization of the total routhian £ w ', 

= £ - A(*|JV|*> - w<*| 4]*), (4) 
leads to the HFB eigenvalue problem for the 2iV x 2iV matrix: 

(5) 

which gives two N x N matrices U and V as eigenvectors and the diagonal matrix E 
of iV eigenvalues which are the quasiparticle routhians (qpr), E^—bijEf '. The complete 
HFB spectrum is composed of pairs of opposite eigenvalues |27| and in Eq. (pj), only one 
eigenvalue from each pair should be selected. The Hartree-Fock hamiltonian h and the 
pairing potential A are constructed from the density matrix p and the pairing tensor k: 





p = V*V T 

K = V*U 



t ' (6) 



according to 



hij — Tij + J^kl Vikjl Plk ~ ~ uJij i ^ 

In order to find a self-consistent solution of the HFB equation @ we have adopted 
a method which relies on the simultaneous use of two different bases of single-particle 
states. The first one is the basis of functions 4>i(r) which converge to the eigenfunctions 
of the Hartree-Fock hamiltonian h and the second one is the canonical basis, i.e., the 
basis of eigenstates of the density matrix p. The iteration procedure is performed in the 
following way. 

Suppose we start at some point with a set of orthogonal functions <f>i(r) and with the 
hamiltonian h in the coordinate representation. First, we perform one imaginary-time-step 
evolution of the set 4>i(r) using h as an evolution operator. The imaginary-time evolution 
|28| has already become a standard technique in the HF coordinate-space calculations. It 
ensures that as iterations proceed, the set <pi(r) converges to the eigenfunctions of h with 
the single-particle routhians ef as eigenvalues, hij—5ij(ef—X). 

Second, we calculate the matrix elements of the pairing potential and the diagonal 
matrix elements hu of the Hartree-Fock hamiltonian in the basis of the new (imaginary- 
time evolved) wave functions <pi{r). As discussed above the matrix converges to the 
diagonal matrix 

and therefore we can safely disregard its off-diagonal matrix elements during the iter- 
ation. On the other hand, for hu^O all the matrix elements Ay are in general nonzero 
both during the iteration and after the convergence is obtained. 
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Third, we solve Eq. (El) to obtain the matrices U and V, the quasiparticle routhians 
Ef as well as the matrices p and k by using Eq. A diagonalization of the density 
matrix, 

Pki = £niW w W5 , (8) 

i 

provides us with the occupation coefficients rii and with the unitary transformation W 
which relates the basis <f>i(r) to the canonical-basis wave functions fi(r)=J2j Wji<f>j(r). 
Finally, we construct the HFB density matrix in the coordinate space 

p(r,r') = Y,n i cp i {r)cp i (r')* , (9) 

i 

which is given as a single sum of products of the wave functions <fi{r). This fact is crucial 
for a numerically effective calculation of local Skyrme densities on the spatial mesh and 
then of various terms in the Hartree-Fock hamiltonian. In fact, at this point the density 
matrix has a form identical to that used in the HF+BCS method although it corresponds 
to an unrestricted HFB variation. At the expense of handling two sets of single-particle 
wave functions we are thus able to solve the HFB equations by using techniques which 
have already been well developed for the HF+BCS method. 

For our purpose, the two-basis method has several advantages. First of all it allows 
us to use the existing techniques to obtain solutions of the cranked HF equations associ- 
ated with a Skyrme interaction directly on a three dimensional mesh. These have been 
presented in refs. p9[ |3(J and have been applied to study the high spin properties of 24 Mg 
and 80 Sr. Because they do not rely on an expansion on a basis, they can describe with 
the same good accuracy the very different shapes that a rotating nucleus can take (some 
examples are shown in Fig. 7 of f30fl). 

Moreover, a separate treatment of the particle-hole and particle-particle channels be- 
comes possible. Although there exist effective forces which have been shown to describe 
reasonably well both channels of interaction simultaneously^ we believe that the present 
knowledge of the pairing effective interactions is not as sophisticated as that concern- 
ing the mean-field. This is especially true as far as its properties relevant to situations of 
time-reversal symmetry breaking such as encountered in high spins physics are concerned. 
Treating in two steps the particle and the quasiparticle spaces, permits us to evaluate the 
gap matrix in the particle basis <f>i(r) with a different interaction than the one used in 
the cranked HF equation. 

Our two-basis method provides also a more natural scheme for limiting the phase 
space than the truncation generated by a fixed basis. Because we are solving the HFB 
problem (|5|) directly in the cranked HF basis, its dimension N can be chosen much smaller 
than that needed for the accurate description of the HF orbitals 4>i(r). Similarly as in 



the phenomenological approaches |37| and based on our former experience, we limit it 
to orbitals for which the particle routhians ef do not exceed a given value e„ iax . This 
value should be much larger than A+A, where A is the physical pairing gap. The results 
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For instance the SkP Skyrme interaction or the earlier and mathematically sounder global 



effective interaction proposed by Gogny |32|, |33|, |34|7 135|, |36 
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presented below have been obtained with e^ ax =A+5 MeV. This leads to a dimension of 
the HFB matrix of the order of hundred, compared to typical dimensions of the cranked 
HF equations of the order of several thousand. 

By applying the two-basis method we have circumvented the difficulties associated 
with the fact that the HFB hamiltonian has a spectrum unbounded from below in the 
coordinate representation. Therefore, a direct application of any variational method (e.g., 
of the imaginary-time method) is impossible. When spherical symmetry is enforced, the 
HFB problem becomes one-dimensional and its solutions can be obtained by explicitly 
solving a differential Schrodinger-like equation. However this route seems to be difficult 
to take for deformed systems which require a three-dimensional treatment. In practice, 
we use the double representation inherent to the two-basis method, to filter those states 
from the unbounded HFB spectrum which are the most physically important. 



2.2 THE SENIORITY PAIRING INTERACTION IN THE CRANKED 
BASIS 

In the present study, we use the simplest possible pairing force: the seniority interaction. 
This force being compatible with the signature symmetry, we can separate the phase 
space in two subspaces of dimensions n\ and n 2 {n\ + n 2 = N) associated with opposite 
signatures. In the following, we denote respectively with latin (i.e. i) and hat-topped 
latin indices (i.e. i) states in the positive and negative signature subspaces. Then the 
non-zero matrix elements of the seniority forces are 

V iJki =-G(j\f\i)(t\f\k) , (10) 

plus all those obtained from the antisymmetry conditions. In fllPf ), T stands for the 
time-reversal operator. This expression corresponds to the lowest-order truncation of the 
pairing matrix elements of a short range force. The pairing strength G is parametrized 
for neutrons (protons) as Gjy = <?at/(H + N) {G z = g z /(ll + Z)), as was done in our 
preceding studies. We discuss later how to select the values of the two parameters and 
9z- 



According to a well established procedure p8| , we consider the matrices U and V in 
the following block form 



u = : " : V = " (11 





in which the first block has dimension n\ and the second dimension n 2 . This leads to 
simplifications for p and k: 



«= ut n ; r = w+ , (12) 






6 



and for h and A: 

where the expressions for h, h and d can be obtained by inserting Eqs. fll"2|) in Eqs.(^). 
Then the HFB equations (pi) reduce to the following iV-dimensional problem: 










(14) 



This last equation specifies our definition of the quasiparticle routhians (qpr) assuming 
that the eigenvectors of the diagonalization problem are ordered according to decreasing 
eigenvalues {n\ first, ri2 last). The first n\ eigenvectors are associated with the positive 
signature. The last ri2 eigenvectors are associated with the negative signature and the 
corresponding eigenvalues are opposite to the negative-signature qpr's. In particular, 
following the discussion in ref. [^] concerning the parity of the number of particles in the 
BCS state, we do not consider the signs of the eigenvalues to decide which qpr's are of a 
given signature. 

In a similar manner, to study one-quasiparticle or two-quasiparticle states, we ex- 
change the U and V components of eigenvectors according to the following procedure. 
First, we choose a given configuration by selecting the label % of the wave function <f>i(r) 
upon which the quasiparticle excitation is to be constructed. Then we scan the ith line 

(Zi V \ 

I looking for the largest component to determine the column (i.e. 

the qpr) for which the exchange of the U and V components must be done. We can also 
determine whether this particular quasiparticle excitation is predominantly of the particle 
or of the hole type; in the former (latter) case the largest component is located in the u or 
u (v or v) blocks. We have found this procedure more stable than that using the values of 
qpr's to select quasiparticle excitations. Indeed, during the course of the HFB iterations, 
the displacement of the Fermi level necessarily associated with the creation of one or two 
quasiparticles makes it difficult to recognize a given configuration solely on the basis of 
the value of its qpr. 



2.3 THE LIPKIN-NOGAMI PRESCRIPTION 

Although the HFB equations contain most of the ingredients essential for the understand- 
ing of the physics of fast rotating nuclei, they present a significant weakness inherent to 
their mean-field nature: for specific values of the angular velocity, they predict phase 
transitions which, in principle, should not occur in a finite system such as the nucleus. 
As we will show, such transitions generate characteristic behaviours in the evolution of 
the moment of inertia versus angular velocity which are incompatible with present data. 
They correspond to a sudden breakdown of neutron or proton pairing correlations and are 
associated with a poor description of particle-number symmetry breaking in situations of 
small particle-number fluctuations. To cure this deficiency while retaining the quality and 
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the simplicity of a mean-field description, the best approach would be a variation after 
projection (VAP) on correct particle number |j27| . Because it appeared too difficult to 



implement it in this first study, we have adopted the approximate projection method pro- 
posed by Lipkin and Nogami p4| , [25| , p6|j . The associated equations which we recall below 



have been tested to provide a good numerical approximation of the VAP in situations 
where both the HFB and the BCS equations predict a collapse of the pairing correlations. 
(See the recent studies of model systems without |?{], ¥1] and with rotation pUfl , and the 
references cited therein.) 

The prescription of Lipkin-Nogami amounts to modify the energy £ (or the total 
routhian £ w in case of cranked HFBLN method) by adding the second-order Kamlah [|43f 
correction: 

£^£-\ 2 (AN 2 ) , (15) 
where (O) denotes the expectation value (\l/|0|\l/). The coefficient A2 is given by: 

(H(AN 2 - (AN 2 ))} - (HAN) (AN 3 ) /(AN 2 ) 

A 2 = 71 —2 ~2 ~2 • ( 16 ) 

(AN ) - (AN ) - (AN) /(AN ) 

For the HFB state, the moments of the operator AA = N — (AN) which appear in the 
definition of A2 are given by 




2tr [ X ] 
4 tr [7x 

3(AA 2 ) +8tr[x(l-6 X )] 




and 



(HAN) = 2 tr [hx] - 3? tr [A«* 7 ] 

(H(AN 2 - (AN 2 ))) =Atr[(h'y + v HF {x})x] - Ktr [Ak*(1 - 8x)] > ( 18 ) 

-tr [u P {7«}(7«)*] 

where the symbol tr [a] stands for the trace of the matrix a and 3ft denotes the real part. 
The mean fields vhf{x} an d vp^n} are the HF potential and the pairing potential 
calculated with the densities x an d 7«, respectively: 

VHFixhj = E Vikjixik ; «p{7«}ij = v m(jK)ki ■ (19) 

A consistent application of the Lipkin-Nogami prescription requires in principle using 
the full effective interaction in Eqs. ([19]). In the present study we calculate A 2 by using only 



the seniority pairing interaction. However, in (18) the effects of the Skyrme interaction 
are still taken into account by the contributions containing the HF hamiltonian h. This 
is consistent with all the other applications of the Lipkin-Nogami prescription performed 
up to date. The fact that particle-number mixing is directly related to the existence of an 
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interaction in the particle-particle channel may be considered as a justification for such a 
procedure. Admittedly, further studies of this question are needed. 

Assuming that the quasiparticle vacuum is a self-consistent solution, one may 
derive another formula for A2 which gives the same value as ([16]) at the end of the iteration 
procedure. In fact, this approach is usually taken in applications, because it gives A2 in 
terms of the two-body interaction only. For the case of the seniority pairing interaction 



such explicit expressions are given in refs. |2Bj and [52| without and with time- reversal 
symmetry breaking, respectively. 

The modification of the HFB equations associated with the Lipkin-Nogami prescrip- 

* 2 

tion is obtained by a restricted variation of A2(AiV ), namely, A2 is not varied although 
its value is calculated self-consistently using fljjp . Based on Eq. (0), this leads [|(J to 
the modification of the HF hamiltonian, h —>■ h — 2A27, while the pairing potential A is 

unchanged. In the present study we have used this method to study 192 Hg, 194 Hg and 

* 2 

194 Pb. However, as noted in ref. [|19j, one may use an equivalent form of (AN )=2tr 
which leads to A — ► A — 4X 2 k and leaves h unchanged. In this way, a better numerical 



stability of the HFBLN equations can be obtained |19[. We have also investigated an 
intermediate prescription: 

/i^/i-A 2 7 ; A ^ A - 2X 2 k , (20) 

which results from inserting the explicit form of the particle- number operator (Nj=5ij) 
only after the variations over p and k are performed. Such a prescription also turns out 
to be numerically stable and was used to study 190 Hg. 



3 Results 

The calculations reported below have been performed using the SkM* parametrization 



45| of the Skyrme force and the seniority pairing interaction as discussed in sec. 2.2. 
The method of solution of the cranked HF equations, which requires breaking of the 
time-reversal symmetry, has been discussed for the Skyrme interaction in refs. [^, |5D| . 
The single-particle wave functions are discretized on a spatial rectangular mesh of points 
spaced by Ax=Ay=Az=lfm. The parity and signature symmetries are assumed to be 
conserved, which allows us to solve the equations by using only points in space which 
have positive values of all three coordinates. 

Although the code does not enforce it, the self-consistent solutions that we have ob- 
tained deviate little from axial symmetry: in all calculations the triaxiality angle 7 never 
exceeds ±0.2°. The average value of the projection of angular momentum on the axis 
perpendicular to the symmetry axis Jq is constrained while that of the quadrupole mo- 
ment Q is unconstrained. In this way, the self-consistent quadrupole moment of the SD 
minimum is found for each value of the angular momentum and the evolution of Q is 
obtained as function of J . 

We use a quadratic constraint on the angular momentum and therefore we are able 
to select values of J along the rotational band equal to even integers. The Lagrange 
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multiplier uj required to obtain a given value of Jo is interpreted as the angular velocity. 
Due to self-consistency the second moment of inertia Sf^ 2 \ defined as 



J* = §, (2!) 



can be expressed in terms of the first derivative of the energy (0) with respect to 



UJ, 



J {2) = ~ (22) 

UJ OUJ 

Neglecting a weak dependence on uj of the Lagrange multipliers A associated with 
proton and neutron numbers, one can also express the angular momentum as the first 
derivative of the total routhian Eq. (f|), 

d£ u 

J ° - -Ir- (23) 

Within the same approximation one therefore obtains an alternate formula for the 
second moment of inertia: 

d 2 S u 
duj 2 

Because formulae (|21|) and involve a first order derivative, they are numerically more 
stable than ([24]) . Moreover, formula (|21 ) which uses the average value of a one-body 



J {2) = —7TT- (24) 



operator is easier to handle than fl2"2"|) which depends on the two-body interaction. We 
have found that an accurate determination of requires a high level of convergence. 
In order to make connection with the experimental method for determining the second 
moment of inertia we have also calculated according to: 

^ Ei+2 — 2£j + £/_2 ' ^ 

where Ej is the total energy £ obtained by constraining the average value of the cranking 
angular momentum to even integers, i.e., Jq=I. 

In practice we have to perform about 200 iterations (or 300 in the crossing regions) with 
an imaginary time-step <5t=0.015fr/MeV. Only then do we obtain the second moments of 
inertia calculated from and (p5|) equal up to 1%. 

At this point one should stress an important difference between the interpretation of 
quasiparticle routhians (qpr's) in the HF and the Nilsson-Strutinsky (NS) approaches. 
In both cases, the value of the angular momentum J is the sum of contributions from 
individual quasiparticle states, Jo=J2iji, where ji are the individual alignments, i.e., 
the average values of the cranking angular momentum J x in the quasiparticle states. 



Therefore, the second moment of inertia (|21|) is also a sum of contributions from individual 
quasiparticle states, equal to the first derivatives of alignments with respect to uj. On the 
other hand, the NS individual alignments are equal to the opposite of the first derivatives 
of the NS qpr's, and hence their contributions to are given by the second derivatives 
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of qpr's. This is no longer true for the HFB qpr's because the HF mean field depends 
self- consistently on the angular velocity (cf. discussion in ref. pOp. Only the total values 
of the angular momentum J and of the second moment of inertia J~( 2 ' are equal to the 
opposite of the first and the second derivatives of the total routhian, Eqs. (|23|) and (|24D, 
respectively. They have, however, no direct relations with the sums of the first and the 
second derivatives of qpr's with respect to u. This fact should be kept in mind when 
analyzing the diagrams of the HFB or HFBLN qpr's. In particular, the changes in 
caused by a quasiparticle excitation cannot not be simply inferred from these diagrams. 

3.1 EFFECT OF DIFFERENT TREATMENTS OF PAIRING 
CORRELATIONS 

Figure [I] presents a comparison of the second moment of inertia calculated within 
the cranked HF, HFB and HFBLN approximations for the nucleus 192 Hg. 

The HF moment of inertiaf] keeps an almost constant value (~ 115fr 2 MeV -1 ) ver- 
sus angular velocity. A similar result obtained within the WSNS method is reported in 
ref. ||22|| . This constancy is in contradiction with the data which show a steady increase 
of from 90 to lSS/^MeV" 1 over the range 100keV</iu;<450keV. One is therefore led 
to conclude that in the mercury region pairing correlations play an important role in the 
superdeformed high-spin states. The increase of J'^ would then be at least partially at- 
tributed to a reduction of the superfluidity as the magnitude of the time-reversal-breaking 
cranking field grows. 

The HFB results obtained with a seniority force = gz=l5. 5MeV do indeed display 
an increase of J 1 ^ (Fig. [I|) which is however too rapid at low rotation frequencies. Despite 
a value of j( 2 \Q)=j( 2 \Uuj=ty which is lower than indicated by experiment, the moment 
of inertia becomes larger than the experimental values in less than ten units of angular 
momentum. In addition, the HFB curve presents two sudden drops at fru;=250keV and 
350keV. These two accidents in the J^ 1 curve are due to the collapses of the proton and 
neutron pairing correlations. This is shown by Fig. |2| which displays the evolution versus 
hu of the proton and neutron pairing contributions to the total energy defined by: 

£?{huj) = -tr [A T <] , (26) 

where the index r refers to neutrons (N) or protons (Z). The fact that the proton pairing 
disappears first and the neutron pairing later is related to the magnitude of the pairing 
correlations at hu=0 which are weaker for protons than for neutrons. This in turn reflect 
the corresponding sizes of the SD shell effects. 

Fig. [| gives the evolution of the neutron and proton single-particle energies versus 
quadrupole deformation in the neighbourhood of the SD well. One sees that for the SkM* 

3 There is a 2fi, 2 MeV _1 reduction of the value of as compared with the calculations of ref. |^0| 
which have been performed with Ax=0.8fm. It measures the absolute inaccuracy of the present calculation 
(performed with Ax=lfm) in replacing the differential operators of the Skyrme functional and of the 
angular momentum by a finite-difference expression on the spatial mesh. 
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interaction, the proton magic gap at Z=80 is large and dominates the neutron shell effect 
at iV=112 which occurs at a slightly larger deformation. WSNS calculations also find 
that this proton gap in the single-particle spectrum is larger than the neutron gap. As a 
consequence, at hu=0 the proton pairing correlations are weaker than neutron ones and 
tend therefore to disappear first. Once the pairing correlations have vanished, the HFB 
moment of inertia becomes identical to that of the HF method. 

The introduction of an approximate treatment of particle number projection by the 
Lipkin-Nogami method modifies the behaviour of and as function of hto. Because 
the projection always leads to an increase of pairing correlations (cf. ref ||46||), we had to 
lower the pairing strengths to gw = gz — 14MeV so as to obtain a moment of inertia close 
to that found with the HFB method at hu=0. Fig. || shows that despite this reduction, the 
proton pairing energy £f (0) is significantly increased over the HFB value while S^(0) is 
much less affected. This reflects another characteristics of the particle-number projection 
which induces larger relative modifications whenever the pairing correlations are weak, 
i.e., when the particle shell effects are large. 

Within the HFBLN method the pairing energies decrease smoothly with the angular 
velocity (Fig. 0). The spurious sudden phase transitions disappear and the calculations 
predict a regular increase of the moment of inertia without any discontinuity (Fig. |l|). Still 
one notes that at tkv—0 the moment of inertia is smaller than the value suggested by an 
extrapolation of the experimental data. Since we consider this work as a first step towards 
a more complete analysis of the effective interaction in the particle-particle channel, we 
have chosen to further renormalize the pairing strength so as to have a value of (0) in 
agreement with experiment. This is obtained by the reduction of both pairing intensities 
to the values of Qn = 9z = 12.6MeV. In all likelihood, such a modification lies within 
the present domain of uncertainty of a pairing force defined entirely on the basis of static 
properties such as quasiparticle energies at ground state or deformation energy curves. 
Dynamical properties such as the moment of inertia will certainly introduce stronger 
constraints. 

With these reduced values of g T , we find that at hu=0 the neutron and proton pairing 
gaps of 192 Hg are equal to Ajv=0.557MeV and A^=0.647MeV, respectively. The evolution 
of pairing energies versus Tluj (Fig. 0) shows an interesting feature: their rate of decrease 
versus Tluj is reduced together with their values. This yields larger pairing correlations at 
high frequencies and hence a slower increase of versus Tluj in better agreement with 
data (Fig. [I]). At the highest spins reached in the present study, the pairing correlations 
are still present for all the nuclei considered. For instance, the pairing energies and 
£f in 194 Pb are at J =52h still of the order of 0.5MeV (see Fig. |). 

Had we been interested in a purely phenomenological analysis of data, we could prob- 
ably have improved the agreement by introducing values of g^ different from those of gz 
while keeping the value of J7"^(0) consistent with the experimental data. On the other 
hand, in our opinion the relative importance of the neutron and proton pairing should 
be ultimately determined by the isospin symmetry. Unfortunately, this symmetry can- 
not be straightforwardly implemented within a seniority interaction which is by defined 
separately for neutrons and protons. We have therefore chosen to defer to further studies 
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such an analysis which will involve a short-range pairing force. All the results, discussed 
hereafter have been obtained with the the pairing strengths Qn = 9z = 12.6MeV. 



3.2 HIGH-SPIN PROPERTIES OF GROUND-STATE SUPER- 
DEFORMED BANDS 

In this section we present results for the ground super deformed bands of the 190 Hg, 192 Hg, 
194 Hg and 194 Pb nuclei. Fig. [| shows the evolution of their charge quadrupole moments 
Q c versus %u. The results for the four nuclei are qualitatively different. The curves 
for 192 Hg and 194 Hg display an expected behaviour (see also Fig. 3 in ref. jl7|): after 
a small increase that can be attributed to the diminishing of pairing correlations, the 
usual effect of quantum rotations [47j takes over and leads to a decrease of Q c . For 194 Pb 



the deformation of the minimum of the SD well at Jo=0 is larger than in the mercury 
isotopes, even when the 4% effect on Q c associated with the two additional neutrons is 
taken into account. This can be explained by considering two characteristics of the Z=82 
subshell effect (see Fig. |3|): first it is shifted to slightly larger deformation than the Z=80 
one; second it is weaker, so that the neutron iV=112 shell effect acts more efficiently and 
pulls the deformation to larger values. Similar properties are visible in the particle proton 



spectrum obtained with the WSNS method (see Fig. 6 of ref. [0). Another specificity of 
the evolution of Q c for 194 Pb is its constant growth versus Tiuj. 

The curve for 190 Hg displays two distinct patterns. At low angular velocities (fru;<325keV) 
the behaviour is similar to that of the other two Hg isotopes. Then occurs a neutron two- 
quasiparticle alignment (see below) so that the neutron pairing is reduced (see Fig. f|). 
This alignment leads to a further reduction of the quadrupole moment as discussed in 



47 



These calculated variations of Q c are smaller than the present level of accuracy of 
the experiments which are still consistent with unchanging quadrupole moments [3^, |49| . 



We hope that coming improvements in the 7-ray detectors will allow the observation 
of such fine differences in the behaviour of the quadrupole moment versus hu. These 
variations are however sufficiently large to significantly modify the dependence of single- 
particle routhians on hu. This is so because deformation effects on single-particle states 
are in general stronger than those induced by rotation. Fig. |6] presents the single-particle 
routhians in 190 Hg calculated along the path of self-consistent quadrupole moments. One 
can see that after the crossing the ftcu-dependence of the routhians is modified, and more 
pronounced because of faster changes of the deformation. 

Fig. |7| shows the single-particle routhians of 192 Hg. There is a large proton shell effect 
for Z = 80 and weaker neutron ones for iV=104, 110, 112 and 118. There exists also a 
weaker proton shell effect for Z=82. This spectrum looks much like that obtained within 



the WSNS |22| , [2311 approach. However, we note some small deviations which may explain 
some of the features that we discuss below. 

In the following, we label the single-particle and quasiparticle routhians by their dom- 
inant components in the asymptotic Nilsson basis. Let us recall that this labelling is 
mostly useful to attribute names to the discussed orbitals. Indeed, the dominant compo- 
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nents have sometimes weights as low as 20% (see ref. ]50| for a detailed analysis). 

In the neutron spectrum, the SkM* force places the intruder orbitals ^[770] 1 ^ 2 , ^[761] 3 y 2 
and i/[752] 5 y 2 about 0.8MeV lower than does the WSNS approach (with respect to the 
rest of the spectrum). On the other hand, the proton spectrum is very similar. The only 
noticeable difference is the existence in our calculation of a bunching of three levels at 
the Fermi surface which leaves the 7r[651] 3 / 2 intruder level isolated with small gaps at 
Z=72 and Z=74. In the WSNS spectrum, these four levels are more evenly spread over 
0.8MeV. The effect of the change of deformation between 190 Hg and 192 Hg can be seen 
by comparing the proton single-particle routhians in Figs. |6] and |7|. At Jo=0 the three 
orbitals below the Fermi level are more spread in 190 Hg than in 192 Hg and the Z=74 gap 
is less visible. This spectrum is therefore more similar to that in ref. \ $% . 

Figure |8| shows the evolution of the second moment of inertia for 192 Hg, 194 Hg 
and 194 Pb. The general behaviour of the data |5j, |7|, [8|, [L3| is correctly reproduced, in 



particular the apparition of a plateau at large htu is obtained. We also reproduce the 
relative positions of the curves for the three nuclei: 194 Hg lies above 192 Hg which itself 
is above 194 Pb. Presently we have no explanation for the slightly too rapid increase of 
the HFBLN moment of inertia. An analogous deficiency is also present in the recent 



WSNS LN calculations |19j for Hg. However, since the results can be very sensitive 



to moderate changes of the pairing force, we think that it could be attributed to a slight 
imbalance between the neutron and proton strengths. The fact that we predict a too 
early occurrence of the plateau for 192 Hg and 194 Hg and the right one for 194 Pb supports 
such a conjecture. 

Our calculations show that the values of J 1 ^ are not simply related to the deformation, 
contrary to what would occure for deformed rigid bodies. For instance, the nucleus 194 Pb 
which, among the studied nuclei, has the largest quadrupole moment, has the smallest 
moment of inertia. This shows to what extent the magnitude of a dynamic quantity such 
as J 1 ^ is unrelated with a static property such as Q c . Our finding for 194 Pb is also in 
agreement with the general observation that for SD nuclei the rate of variation versus 
Tiu) is correlated with the occupation of the intruder orbital 7r[642] 5 / 2 . 

The calculated and experimental |23| second moments of inertia of 190 Hg (Fig. ||) grow 
also slowly as a function of the angular momentum until fou;=325keV where the calculation 
gives a peak due to the two-quasiparticle alignment associated with the ^[761] 3 / 2 orbital, 
see Fig. Note that it is precisely at this frequency that the single-particle routhianfj] 
^[761]3 //2 crosses the Fermi level (Fig. P). This illustrates the fact that in the HFB theory 
the occupation probabilities are not directly determined by the eigenvalues of the mean- 
field hamiltonian. Although the experimental error is still large, the last point of the data 
suggests that a crossing will happen at a slightly higher angular velocity. 



4 The superscripts +, — or ± at the Nilsson labels denote the signatures of orbitals 
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3.3 QUASIPARTICLE ROUTHIANS IN THE HG REGION 

In the nuclei which differ from 192 Hg by no more than two nucleon units, 30 SD bands 
have already been observed (in 191 Au, 190 " 4 Hg, 192 " 4 T1 and 194 Pb g, |, § g, @, ft |, |, 0, 
11, 12, |l^] From this information some general features have already become apparent .0 

• It appears that pairing correlations are always playing a significant role. Except for 
two bands in 192 T1, the of all considered SD bands is growing. 

• The nuclei having a comparable content of intruder states have similar slopes of 

as function of hu. The fact that this is independent of whether the nuclei are odd or 
even is an additional indication that dynamic pairing correlations (as described for 
instance by particle number projection) are present and are at least as important 
as the static ones contained in a simple HFB description. 

• As in the rare earth region, there are also specific trends associated with the content 
of the proton or neutron intruder-state configuration. For instance, one can group 
several SD bands into a "family" associated with the ground SD band in 192 Hg. 
These bands have the fastest growing second moment of inertia . Among them 
there are several examples of pairs of identical bands. 



Recently another family associated with two SD bands in 193 T1 has emerged [IT]. It 
includes mostly bands assumed to have the 7r[642] 5/ , 2 orbital blocked. One can also 
add to this family the ground SD band in 191 Hg which is interpreted as having the 
z/[761] 3 / 2 orbital blocked. If this orbital assignments are correct, this would indicate 
that the effects of proton and neutron intruder states are of a similar magnitude. 

A few bands, one recently discovered in 193 Hg |9[] and two in 192 T1 |2|, [TI], display 
a markedly different behaviour. The latter two can be explained by a blocking of 
two (neutron and proton) intruder states at the same time. Similar bands have also 
been observed in 195 T1 ||11|| . 



We first consider the HFBLN quasiparticle routhians for 192 Hg (Figs. [11] and |T2"p . 
Because we are going to discuss results for neighbouring nuclei, we have labeled the 
routhians (at hu=0) not only with their dominant Nilsson configurations but also with a 
letter p or h to indicate whether they are built on a particle or a hole configuration. If 
a nucleus is related to 192 Hg by an addition of particles, the Fermi level will raise. One 
expects that the qpr's built on hole states will be pushed up while those built on particle 
states will go down. The opposite effect will occur when the nucleus is obtained by a 
subtraction of nucleons. 



5 Because this is a rapidly evolving domain of physics, there is unfortunately no consensus on the 
labeling of the bands in the available literature. Sometimes signature partners are denoted with consec- 
utive numbers, sometimes with the same number and alphabetical (a and b) indices. We shall therefore 
avoid such labeling and when necessary, refer to the bands by what seems to be their most probable 
structure according to today's models. 
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In the proton diagram (Fig. |TT]), the qpr associated with the 7r[642] 5//2 orbital is well 
separated from the other ones. It displays a rather gentle alignment pattern and a very 
small signature splitting. This is in agreement with the observation that in all Tl isotopes, 
the bands can be sorted in pairs of signature partners. Moreover, the fact that their 
moment of inertia is always growing less rapidly than that of 192 Hg is consistent with 
having one intruder orbital blocked. 

Experimental routhians can be extracted from the information for the two bands of 
193 T1 by taking 192 Hg as a reference and making reasonable assumptions for the spins. 
They show a decrease of this particular qpr of about 0.4MeV in the range of angular 
velocities of 110<7iu;<360keV and a maximal signature splitting of 30keV. In our cal- 
culations we find 0.6MeV and 80keV, respectively, for these two quantities. A similar 
analysis of the 7r[642] 5 / 2 qpr in 194 Pb leads to smaller numbers: 0.52MeV and 50keV. It 
therefore seems that our calculations predict a slightly too large proton alignment. This 
could explain why the HFBLN moments of inertia (Fig. |3|) increase somewhat faster than 
the experimental data. 

The qpr's associated with the three hole levels 7r[532] 3 y 2 , 7r[530] 1 / 2 and 7r[411]-,y 2 are 
located around lMeV. For lighter neighbouring elements, these qpr's become the lowest 
ones and replace in the qpr spectrum those built on the particle orbital 7r[642] 5 ^ 2 . In 
191 Au our calculation would therefore predict three bands close in energy. The one built 
on 7r[411] 1 ^ 2 would not have an easily detectable partner because of the strong signature 
splitting. It is tempting to say that this is the band which has already been observed. 
The other two result from a strong interaction between the levels 7i [532] 3 y 2 and 7r[530]-,y 2 , 
which pushes down the qpr's 7r[532]3 / , 2 and 7r[530]w 2 of opposite signatures. 

In this respect, our calculation differs slightly from that using the WSNS method. 
Indeed, because the particle routhian built on the orbital 7r[532] 3 ^ 2 is lower in the spectrum 
of ref. H], there is less interaction with the one built on 7r [530] 1 y 2 . Since this orbital has a 
large signature splitting, one should then in principle expect to find only one other band 
without a partner. However, we note that the signature splitting of levels (see for instance 
^[411)^2 in Fig. 3 of ref. [|J) in the WSNS approach is smaller than the one found in the 
present work. Therefore their predicted signature partners are never very far in energy. 
In addition, the negative-signature intruder orbital 7r[651]3^ 2 obtained in ref. seems 
to interact earlier (fru;^250keV) with 7r[411]^ 2 . It crosses it and appears at the Fermi 
surface at /ic<j=0.300keV. In our calculation, until ft,co>=410keV, the 7r[651]3^ 2 state is still 
not favoured as compared to the orbitals 7r[411]^ 2 and 7r[532] 3 t J /2 . Further experiments on 
gold isotopes may certainly help clarify this situation. 

Let us now consider the qpr spectrum for neutrons (Fig. [12|). At Jo=0 one first finds 
a group of three levels. Two are of the particle type, i/[512] s / 2 and i/[624] 9 / 2 , and one of 
the hole type i/[642] 3 / 2 . Further up there is a group of four levels. Among them there 
are two routhians built on intruder orbitals: the particle- type z/[752] 5 ^ 2 and the hole-type 
z/[761] 3 / 2 . Because of a strong interaction between these two intruder routhians, they 
almost immediately exchange their character when hu increases. In particular, it seems 
more appropriate to assign the steep downgoing branch which comes out of the qpr labeled 
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z/[752] 5 / 2 to the orbital ^[761] 3/ , 2 . At this point, we will only justify this affirmation by 
noting that the signature splitting of the latter routhian is larger than that of the former 
one (see Fig. |7|) because it has a smaller value of K. Later, we will analyze this point 
in more details when comparing the neutron qpr spectra of the four mercury isotopes 
studied in this work. 

If we try to extrapolate from this spectrum that of 191 Hg, we expect that the most in- 
teresting configurations will be those associated with the hole-type qpr's: those of the two 
signature-partner states of 2/[642] 3 , 2 and the negative-signature down-going state z/[761] 3 , 2 . 
We note that for /iu;>300keV the latter qpr is the least excited one. Moreover, in 190 Hg the 
same state is the least excited already at /k<j=120keV (Fig. |10D . One can therefore infer 
that the corresponding crossing will happen in 191 Hg at an intermediate angular velocity. 
Since the SD bands are essentially fed at high spins, it is natural to expect the ^[761] 3 / 2 
band to be the most populated. Then should come two bands with moderate signature 
splitting associated with z/[642];y 2 . This corresponds well to the experimental situation. 
The extraction of the experimental qpr's for the two signature-partner bands shows an 
approximate overall decrease of 150keV with a maximal splitting of about 40keV over the 
range of 150</ta;<380keV. For these two numbers we find 330keV and lOOkeV, respec- 
tively. Therefore, as was the case for protons in 193 T1, we find also a neutron alignment 
which is slightly too large. 

We consider now the heavier isotope 193 Hg. A hole state such as z/[642] 3 / 2 will now 
become higher than the two particle states z/[512] 5/ , 2 and i/[624] 9//2 . On the latter orbital 
two SD bands with a negligible signature splitting can be built. Candidates for these 
bands exist in the experimental data. Similarly, there should also be a band built on 
the z/[512] 5/ 2 qpr. This band should be very similar to the two previous ones because of 
the strong similarity of their qpr's which remain almost parallel to each other. This has 
also been observed. In addition our diagram shows that the negative-signature orbital 
i/[512]^ 2 interacts strongly with zy[761]^y 2 . Such a phenomenon has precisely been seen 
in the experiment [|51], [|, []]. Moreover, recently a candidate for the i4761] 3 / 2 band has 
been observed with a reduced J 1 ^ which is compatible with an upbending behaviour 
of the associated qpr, as is visible in our qpr spectrum. The Tiuj value for the crossing 
is measured to be equal to 250keV and compares favourably with the value of 200keV 
obtained in our calculations. 

Let us now pursue this discussion by considering in more details the evolution of the 



neutron qpr spectrum as function of the neutron number for the mercury isotopes. Fig. 13 



includes selections of levels from the qpr spectra obtained for the ground SD bands of 
190 Hg, 192 Hg and 194 Hg. We have also included the spectrum obtained by a calculation 
of 193 Hg treated as an even nucleus, using the fact that the number of particle of a BCS 
wave-function is only defined on the average. The qpr's shown in Fig. [13] for 193 Hg result 
from a calculation in which we have arbitrarily imposed that the average neutron number 
is equal to 113. This allows us to study the influence of the position of the Fermi level on 
the crossing of the negative-signature ^[512] 5 / 2 states with the N=7 intruder states. 
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The qpr's shown in Fig. [13] are z/[642] 3 / 2 and ^[761] 3 / 2 for the holes states, and z/[512] 5/ , 2 
and i/[752] 5 y 2 for the particle states. To make the figure more readable, we have eliminated 
the qpr's associated with the hole states ^[505] 11//2 and ^[640] ^ which play no role in 
our discussion and the particle state z/[624] 9 / 2 . The latter qpr shows no sign of signature 
splitting. In all Hg isotopes, it remains parallel to the ^[512]^ 2 orbital with an energy shift 
which decreases gradually from 150keV in 190 Hg to almost zero in 194 Hg. As mentioned 
above, two signature-partner bands in 193 Hg have been assigned to this particular orbital 
i/[624] 9/2 . 

The 190 Hg qpr diagram shows that the single-particle dynamics is determined mostly 
by the alignment of the strongly signature-split orbital z/[761] 3//2 whose negative-signature 
partner crosses the t/[642] 3 / 2 & t hu=120keV , then becomes negative at hu=300keV and 
generates an alignment at hcu=325keV . At this frequency, the experimental indicates 
no alignment but the last available data point suggests a possible crossing at a higher 
angular frequency. Moreover, we obtain a slightly too large slope of the second moment 
of inertia prior to the peak. Both facts suggest again that our calculations predicts an 
alignment faster than seen in the experiment. 

From the separation of the qpr's, we estimate that the interaction between the ground 
SD band and the two-quasiparticle band is approximately 300keV. According to our di- 
agram it seems also that a negative-parity band built on the z/[642] 3 / 2 will become yrast 
beyond frcu=400keV. 

In the qpr diagrams of 192 Hg, 193 Hg and 194 Hg (Fig. [13]), there is a modulation of the 
interaction between the intruders ^[761] 3 ^ 2 and z/[752]^ 2 , and the z/[512]7, 2 orbital as the 
hole- type qpr's move up in the spectrum. In 194 Hg, the ^[761] 3 / 2 is pushed so high by 
the displacement of the neutron Fermi level that despite its rapid descent, its negative- 
signature branch z/[761] 3 / 2 does not interact with u[752]^ 2 before huj=3Q0keV . Had we 
pursued the calculation to higher spins, an extrapolation suggests that an interaction with 
the i/ [5 12] would have occurred near hu=440keV . It is interesting to see that in these 
four isotopes of mercury, the i>[752]^ 2 curve displays precisely the upbending behaviour 
of the experimental routhian of the band recently discovered in 193 Hg ||. Qualitatively 
speaking, one could then say that the SD band interaction observed in 193 Hg concerns 
the i/[761]^ 2 and i/[512]g/ 2 states while the last SD band in the same nucleus is built on 
the z/[752]^, 2 state. Of course, we must keep in mind that within the energy range of the 
order of the pairing gap around the Fermi level, one cannot distinguish between holes and 
particles. As shown in the qpr spectrum of 190 Hg (Fig. |TJj) the hu=0 identity of the qpr's 
is rapidly lost if they happen to be close in energy and have large Coriolis coupling matrix 
elements, as is the case for z/[752] 5 / 2 and i/[761] 3 / 2 . 

The qpr spectrum of 193 Hg is intermediate between that of 192 Hg and 194 Hg. One 
sees that the lowering of the ^[761] 3//2 is not large enough to induce an interaction at the 
values of hu where it has been observed. In our calculation the interaction only occurs 
at hu=380keV . Nevertheless, there does not seem to be any need to invoke dynamical 
octupole correlations to explain the value of hu at which the band crossing occurs in 
193 Hg p2| |. This is also in line with the results of the Generator Coordinate Calculations 



18 



53J which find that a collective one-octupole-phonon SD band would be excited by 
approximately 1.7MeV with respect to the ground SD band, and therefore difficult to 
detect with presently available techniques. In fact, a mere 200keV up shift of the Jo=0 
position of the ^[761] 3 / 2 intruder state would certainly eliminate most of the discrepancies 
between our calculations and experiment. We will come back to this question in the 
concluding section. 

The neutron qpr spectrum of 194 Pb is presented in Fig. [TJ]. A comparison with the 
isotone nucleus 192 Hg (Fig. [12]) shows the influence of the proton-induced polarization. 
According to the above discussion, because the Z=82 SD shell effect is weaker than the 
Z=80 one, the deformation of 194 Pb is larger. As a result (see for instance Fig. |3j), the 
z/[752] 5 / 2 orbital is approximately equidistant to z/[761] 3 / 2 and z/[512] 5 / 2 at hu=0 (in 192 Hg, 
it is closer to ^[761] 3//2 ). Now instead of two successive interactions, one observes a global 
interaction involving the three negative-signature levels for values of Tiuo ranging from 
150keV to 450keV. This behaviour so different from that observed in 192 Hg results from 
a mere 50keV shift of the i/[752] 5 , 2 qpr away from the ^[761] 3 / 2 orbital. 



4 Summary and conclusions 

In order to analyze nuclear structure at high angular momentum, we have constructed a 
cranked Hartree-Fock-Bogoliubov code. We solve the equations iteratively by performing 
in sequence the following operations. In the first one we evolve the single-particle wave 
functions by one imaginary-time step using the cranked Hartree-Fock mean field. In the 
second one we build the HFB matrix and solve the associated eigenvalue problem. After 
selecting the occupied quasiparticle states we compute the one-body density matrix and 
the pairing tensor. Then the eigenf unctions of the density matrix and the associated 
occupation probabilities are used to define the potential and the effective mass which 
enter the cranked Hartree-Fock mean field. This completes one iteration. 

This method of solution which requires the simultaneous handling of two single-particle 
bases, allows an easier separation of the particle-hole and particle-particle (pairing) chan- 
nels. It also permits a significant reduction of the size of the HFB problem. Since it is our 
feeling that the mean-field part of the effective force is reasonably well known (at least 
in the vicinity of the stability line), we plan to use this method to explore in more detail 
the pairing channel whose influence on the rotational properties of nuclei is expected to 
be important. 

A HFB test calculation for the nucleus 192 Hg has shown the limitations of the mean- 
field method. At the angular velocities for which either proton or neutron pairing disap- 
pears , the dynamical moment of inertia displays abrupt decreases which are not compat- 
ible with the experiment. It is known that this deficiency results from a poor treatment 
of pairing correlations by the mean-field approach in the critical and subcritical regimes. 
It is also known that this can be improved by a variation after particle number projection 
(VAP). We have used the Lipkin-Nogami prescription (HFBLN) which provides a good 
approximation of the exact VAP. As expected, the HFBLN results for do not display 
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spurious phase transitions and are in good qualitative agreement with data. 

We have studied the ground super deformed bands of the four nuclei 190 Hg, 192 Hg, 
194 Hg and 194 Pb. The agreement is generally good. Using a single parameter for the 
strength of the pairing force, and adjusting it to reproduce the experimental dynamical 
moment of inertia at low frequency, we can obtain the slow upgoing behaviour of data 
including the appearance of a plateau around hu=300keV. We also find that the of 
194 Pb is smaller than that of 192 Hg which itself is smaller than that of 194 Hg. An analysis 
of the quasiparticle routhians has shown that our spectrum is compatible with the single- 
particle structure of the SD well as it can be deduced from the present experimental 
knowledge based on 30 bands in nuclei which are close neighbours of 192 Hg. This includes 
the interesting band-interaction effect observed in 193 Hg. 

Despite this good agreement, we have detected some weaknesses in our results. The 
most evident one is a slightly too fast alignment. It seems reasonable to assign this 
property to a too simple form of the pairing interaction. Several lines of improvements 
appear possible. First one could consider modifying the relative balance of the neutron and 
proton strengths. One can also study the effects coming from a more refined interaction 
such as for instance a density-dependent zero-range force. 

It has also become apparent that a good description of data requires a very precise 
determination of the particle spectrum at tkv—0. We have seen that our spectrum differs 
from that obtained within the WSNS approach ]17| as far as the positions of some neutron 
intruder states are concerned. The differences between the models may originate from 
small variations on the position of the intruder states at sphericity. Indeed any such 
difference will be emphasized by the large deformation associated with superdeformation. 
Another cause might be the value of the equilibrium deformation at the minimum of the 
SD well. Since we are often dealing with states which in the Nilsson diagram are rapidly 
up- or downgoing, a modification of the equilibrium deformation (of the order of le.b for 



Q c ) could explain the differences between our spectrum and that of ref. [|17 |. Since we 
have seen that the deformation of proton and neutron shell effects are slightly offset, a 
small modification of proton versus neutron pairing strengths could probably realize this 
effect. In fact, the uncertainty as to the exact magnitude of the deformation is only one 
more example stressing the importance of polarization effects on the high spin behaviour 
of nuclei: in this work, we have also discussed a case where a change in proton occupation 
modifies the behaviour of the neutron quasiparticle routhians. 

Altogether the differences between various approaches in predicting the high-spin 
structure remain rather small, at least for those which also successfully describe col- 
lective properties at zero spin. Such differences could even be considered as being within 
the predictability range of the methods themselves. Nevertheless the challenge is there: 
the quality of present experimental information has introduced us into an era where much 
stricter requirements will be imposed on nuclear structure models. 
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Figure Captions 



Figure 1: The second moment of inertia of the ground SD band of 192 Hg. The data 
(open circles) are taken from ref. ||. The other curves correspond to the HF (solid line), 
HFB with # r =15.5MeV (triangles), HFBLN with ^ T =12.6MeV (solid circles) and HFBLN 
with g T =14MeV (solid squares) calculations. 



Figure 2: Pairing contribution to the total energy of the ground SD band of 192 Hg. The 
curves correspond to the HFB with g T =15.5MeV (dashes), HFBLN with g r =12.6MeV 
(solid line) and HFBLN with g T =14.MeV (dots) calculations. 



Figure 3: Nilsson diagram obtained by a quadrupole constrained HFBLN calculation for 
the nucleus 192 Hg at hu=0. The deformation is measured by the charge quadrupole mo- 
ment in e-b. The levels are indexed by their dominant asymptotic configuration. Positive- 
and negative-parity levels are drawn as solid and dot-dashed lines, respectively. The dot- 
ted line indicates the Fermi level. The arrows indicate the charge quadrupole moment of 
the SD minimum. 



Figure 4: Proton (dashes) and neutron (solid line) pairing energies in MeV versus angular 
velocity for the ground SD bands. 



Figure 5: Charge quadrupole moments in e-b versus angular velocity for the ground SD 
bands. 



Figure 6: Single-particle routhians as function of the angular velocity for the ground SD 
band of the nucleus 190 Hg. The convention for the different (parity, signature) combi- 
nations is: (+, +) solid line, (+, — ) dashed line, (— , +) dot-dashed line, (— , — ) dotted 
line. 
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Figure 7: Same as in Fig. ^ for 192 Hg. 



Figure 8: The second moments of inertia J 1 ^ of the ground SD bands of 194 Hg, 192 Hg and 
194 Pb. Data for 194 Hg P2"| , 192 Hg [|] and 194 Pb [0] are indicated respectively by empty 
squares, circles and triangles. HFBLN results are drawn as dashed ( 194 Hg), solid ( 192 Hg) 
and dot-dashed ( 194 Pb) curves and solid symbols. 



Figure 9: The second moment of inertia of the ground SD band of the nucleus 190 Hg. 



Data [23[] are indicated by open squares. HFBLN results are drawn as a solid curve. 



Figure 10: Neutron quasiparticle routhians as function of the angular velocity for the 
ground SD band of the nucleus 190 Hg. The convention for the different (parity, signature) 
combinations is: (+, +) solid line, (+, — ) dashed line, (— ,+) dot-dashed line, (— , — ) 
dotted line. 



Figure 11: Same as in Fig. 10 for protons in Hg 



Figure 12: Same as in Fig. [10] for Hg. 



Figure 13: Comparison of the evolution of selected neutron quasiparticle routhians as 
function of the angular velocity for the ground SD bands of the nuclei 190 Hg, 192 Hg, 193 Hg 
(see text) and 194 Hg. The convention for the different (parity, signature) combinations is 



the same as in Fig. 10 



Figure 14: Same as in Fig. |10| for 194 Pb. 
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